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Abstract 

A  numerical  analysis  Is  made  of  developing 
and  developed  turbulent  flow  In  the  annular 
region  between  a  cylinder  moving  at  constant 
velocity  within  a  fixed  concentric  tube. 

Turbulent  shear  Is  modeled  by  eddy  v 1 scosl ty .and  a 
uniform  velocity  Is  assumed  at  the  entrance  to  the 
annular  region.  The  computations  extend  and 
modify  the  method  of  Sud  and  Chaddock  (1981)  to 
arbitrary  Reynolds  numbers  and  radius  ratios. 
Approximate  formulas  are  given  for  both  developing 
and  developed  values  of  pressure  drop,  shear 
stress  on  the  inner  and  outer  walls,  and  total 
cylinder  drag.  _ 

Nomenclature 


C  •  parameter  defined  in  Eq.  (11) 

Dh  -  hydraulic  dlamater,  2(r0-ri) 
f  -  friction  factor,  Eq.  (’0) 

F  -  drag  force,  Eq.  (19) 

L»  -  entrance  length,  Fig.  7 

n  *  eddy  viscosity  parameter.  Table  1 

p  -  pressure 

r  -  radial  coordinate 

rmt  "  radius  of  maximum  velocity,  Fig.  2 

Rev  •  Uvri/v 

Reb  -  UpDh/y 

u  *  axial  velocity 

lip  *  average  annular  velocity,  Eq.  (1) 

Uv  «  inner  cylinder  or  vehicle  velocity 

U  -  lnvlscld  core  velocity,  Fig.  2 

v*  -  friction  velocity,  (ty/p)1'2 

x  *  axial  coordinate 

y  -  coordinate  normal  to  wall,  Fig. 2 

y*t  -  eddy  viscosity  crossover  point.  Table  1 

4  -  boundary  layer  thickness 

p  -  density 

v  -  kinematic  viscosity 

t  -  shear  stress 

T»  *  2t/pUv2 
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law-of-the-wall  coordinates,  Eq.  (9) 


Subscripts 


inner  wall 
outer  wall 
either  wall 
fully  developed 
lnvlscld  core  value 
In  the  entrance  region 
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The  problem  of  turbulent  flow  in  a  AoMKMfri'J  ■ 
annulus  has  been  of  Interest  for  at  least  eight 
decades,  beginning  with  the  experimental  results 
of  Becker  (1907).  Numerous  papers  have  treated 
the  problem  of  pressure-driven  flow  between  fixed 
cylinders.  For  fully  developed  flow,  the  eddy 
viscosity  theory  of  Quarmby ( 1968 )  is  In  good 
agreement  with  his  experimental  results  (Quarmby 
1966).  For  both  laminar  and  turbulent  flow  in  a 
fixed  annulus,  the  point  of  maxi  ,  axial  velocity 
rmt  Is  closer  to  the  Inner  radius  (Lawn  and  Elliot 
1972). 

Developing  annular  flow  is  less  well 
documented,  but  an  integral  analysis  by  Wilson  and 
Medwell  (1971)  shows  a  development  length  for 
turbulent  flow  of  about  ten  hydraulic  diameters 
for  Reynolds  numbers  from  1011  to  3x10^  and  radius 
ratios  from  1.25  to  5.0.  Their  results  agree 
reasonably  well  wir'.  data  by  Okllshl  and  Serovy 
(1967). 

The  problem  under  study  here  is  the  annular 
flow  driven  by  uniform  axial  motion  of  the  tnner 
cylinder,  which  simulates  the  motion  caused  by 
vehicles  travelling  In  tubes  or  by  the  launching 
of  torpedoes  or  cylindrical  buoys.  Experimental 
data  on  vehicles  in  tubes  are  summarized  by 
Davidson  (1979),  and  theoretical  results  are  given 
in  the  book  by  Hammlt  (1973).  Implicit  In  the 
analysis  of  the  annular  flow  induced  by  a  moving 
cylinder  Is  the  assumption  that  the  volume  swept 
out  by  the  vehicle  passes  through  the  annulus  to 
fill  the  void  behind  the  vehicle,  as  in  Figure  1. 
This  simulates  a  very  long  outer  tube,  whereas  In 
a  short  tube  the  inlet  and  outlet  can  serve  as  a 
source  and  sink  of  fluid,  with  less  volume  passing 
through  the  annulus.  The  only  detailed  analysis 
of  the  flow  shown  In  Figure  1  is  given  by  Sud  and 
Chaddock  (1981),  for  both  developing  and  developed 
flows.  After  presenting  a  general  analytical 
technique,  Sud  and  Chaddock  give  computed  results 
for  only  a  single  case,  r0  ■  8  ft  (2.99  m),  r(  *  6 
ft  (1.83  m),  and  Uv  -  200  ft/s  (61  m/s),  at  two 
air  densities  (Reynolds  numbers). 

It  Is  the  purpose  of  the  present  paper  to 
extend  the  analysis  of  Sud  and  Chaddock  to 

arbitrary  radius  ratios  and  Reynolds  numbers,  for 
both  developing  and  developed  flow.  The  developed 
flow  analysis  Is  identical  to  Sud  and  Chaddock. 

For  developing  flow,  corrections  are  made  to  their 
analysis,  as  it  Is  felt  that  Sud  and  Chaddock’s 
results  in  the  entrance  region  are  somewhat 
inaccurate. 
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Figure  1.  If  a  vehicle  moves  through  a 

long  tube,  the  volume  swept  out 
passes  through  the  annulus. 


Fully  Developed  Flow 

The  geometry  of  tne  flow  ia  shown  In  Figure  2. 
With  coordinates  fixed  to  the  vehicle,  both  the 
upstream  fluid  and  the  outer  t  :be  wall  appear  to 
move  at  vehicle  velocity  Uv.  Based  on  the 
displacement  concept  of  Figure  1,  the 
incompressible  continuity  relation  requires  that 
the  average  velocity  Up  In  the  annular  region  be 
related  to  vehicle  velocity  by 


♦  3  ♦  ♦  _  2  *  * 

dy  l*nuy[l*  exp(-n  u  y  )] 

The  turbulent  core,  y  ^  <_  y  ^  yml,  ,s  modeled  with 
Karman's  similarity  hypothesis  for  eddy  viscosity: 


?  ♦  *  »  2 
d~u  <  (du  /dy  ) _ 

.2  *  ♦  1/2 

dy  ( t/tu  -  du  /dy  ) 


<«  0.36 
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In  both  cases,  the  local  shear  stress  is  related 
to  local  radius  through  a  momentum  balance: 


,2  2  , 

r  (r  -r  ) 
w _ mt 

,  2  2  , 

r  ( r  -  r  ) 
w  mt 


(7) 


At  high  Reynolds  numbers,  Deissler’s  damping 
parameter  n2  equals  0.0151!  with  the  cross-over 
point  y*j  *  15.0.  At  lower  Reynolds  numbers,  values 
of  n2  and  were  taken  from  a  graph  given  by 
Quarmby  (19b8).  Some  typical  values  are  given 
here  in  Table  1.  These  values  are  not  thought  to 
be  especially  accurate  but  are  retained  for 
comparison  to  the  analyses  of  Quarmby  (1968)  and 
Sud  and  Chaddock  (1981). 


Equations  (5,6,7)  are  a  set  of  differential 
equations  to  solve  for  the  velocity  profiles  u( 
and  u0  in  the  Inner  and  outer  regions.  From  Fig. 
2,  the  boundary  conditions  are 


Up  -  Uvb2/(b2-l)  (1)  . 

where  b  -  r0/r(  is  the  principal  geometric 
parameter  in  the  analysis. 

In  the  fully  developed  region,  u  -  u(r)  only, 
and  the  flow  may  be  divided  Into  an  Inner  region 
"1"  and  an  outer  region  "o" ,  separated  by  the 
point  of  maximum  velocity,  rrat.  Since  there  is  no 
acceleration  in  the  developed  region,  a  momentum 
balance  shows  that  the  inner  and  outer  wall  shear 
stresses  are  related  by 


Ti  b(a2-l) 
where  a  *  rmc/ri  . 


(2) 


There  are  two  characteristic  Reynolds 
numbers:  for  the  vehicle  motion,  Rey  -  Uvri/v, 
and  for  the  annular  flow,  Rep  -  UpDp/v,  where  Dp  - 
2(r0-rt)  is  the  hydraulic  diameter.  The  two  are 
related  by  geometry  and  continuity: 

Reb  -  2b2Rev/(b*l)  (3) 


The  analysis  follows  Sud  and  Chaddock  (1981). 

Both  the  inner  and  outer  regions  are  divided  into 
two  computational  regions  as  in  Fig.  2:  (1)  the 
sublayer,  and  (2)  the  turbulent  core.  The  flow  Is 
assumed  steady,  incompressible,  and  turbulent, 
with  impermeable  walls.  The  variables  are 
non-dlmenslonallsed  in  law-of-the-wall  fashion: 

» 

u*  *  u/v»  :  y*  -  yv»/v  ,  v*  -  (tu/p)1/2  (**) 


The  sublayer,  0  <  y  <  y*.  is  modeled  with 
Delssler's  eddy  viscosity  formulation: 


yt 

*  9.  : 

'■>1 

-  0 

(8) 

y0 

•  0  ; 

,Jo 

■  Uv 

The 

equations 

are 

solved  by  an  iterative  technique 

similar  to  Sud  and  Chaddock  but  modified  by  Koclow 
(1985).  One  begins  by  guessing  tq  and  a  (or 
rmt),  whence  i[  can  be  computed  and  the  velocity 
profiles  uj  and  u0  computed  by  numerical 
(Runge-Kutta)  Integration  of  Eqs.  (5,6,7).  The 
value  of  t0  is  gradually  adjusted  and  the 
integrations  repeated  until  the  core  velocity 
Uj0  exceeds  Up.  If  the  two  profiles  do  not  meet 
at  rmt  with  Usi  *USo  then  a  new  value  of  a 
is  computed  from  the  relation 


a  .  (b(I  *  a  -  u;e/U;t  (9) 

b  ♦  a 
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Figure  2.  Velocity  and  boundary  layer 
development  of  fluid  in 
the  annular  region  (modified 
from  Sud  and  Chaddock  1981). 
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The  entire  process  ts  tner.  repeated  until  the 
value  of  a  converges.  The  profiles  are  then 
integrated  to  determine  the  average  velocity  Uj-, , 
which  1 s  compared  with  Eq.  (1).  If  there  Is  a 
mismatch,  the  value  of  t0  is  again  adjusted  and 
the  entire  double-iteration  process  is  repeated 
until  final  convergence.  With  the  profiles  known, 
the  pressure  drop,  friction  factor,  and  other 
parameters  are  computed. 


Table  1  Eddy  Viscosity  Parameters  at  Low  Reynolds 
Numbers 
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Using  the  above  procedure,  Kotlow  (1985)  ha3 
tabulated  fully-developed  flow  parameters  for  a 
range  of  Rev  •  lO^-lO®  and  radius  ratios  b  » 
1.01-2.0.  The  fully-developed  friction  factor  Is 
deftned  by 
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fd 


(-dp/dx)  D. 
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8(r  t  ♦  r 
o  o 

(ro*r1)  pu 
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(10) 


Values  of  ffj  are  plotted  In  Fig.  3  from  the  above 
theory.  The  pressure  drop  is  very  large  at  small 
values  of  b  due  to  the  high  bulk  velocity  In  the 
annulus,  Eq.  (1).  The  curves  In  Fig.  3  can  be 
collapsed  approximately  Into  a  Moody-chart  type  of 
curve  by  scaling  Reynolds  number  and  pressure  drop 
with  U(j  rather  than  Uv.  Note  the  agreement  in 
Fig.  3  with  the  two  computed  points  of  Sud  and 
Chaddock  (1981)  for  b  -  1.33. 


The  computed  curves  In  Fig.  3  may  be  fit  to 
the  following  algebraic  formula  with  good 
accuracy: 


Figure  3.  Theoretical  fully-developed 
friction  factors  in  the 
annular  region. 


In  the  fully-developed  region,  the  drag 
experienced  by  the  vehicle  would  be  a  combination 
of  pressure  difference  and  inner  wall  shear.  For 
a  cylinder  of  length  L,  the  drag  would  be 


Ffd  -  ,r1L(  2t,  -  r,  £  > 


(IV) 


This  estimate  must  be  supplemented  by  values  of 
nose  drag,  wake  drag,  and  the  Increased  pressure 
drop  and  shear  In  the  developing  region  (Davidson 
1979). 


* 


ffd  "  c?  [0.001  ♦  2.8/(log10Re*)3-l],  (11) 


where  C  -  — b  —  and  Re*  -  Re„ 

2(b2-l)  b+1  V 

The  error  when  compared  to  Fig.  3  Is  about  ±3*  in 
the  range  Rev  -  1(P-10®. 


Finally,  the  moving  Inner  cylinder  causes 
the  point  of  maximum  velocity  to  shift  toward  the 
outer  wall,  unlike  the  fixed  annulus  case.  This 
position  Is  shown  in  Fig.  9  and  becomes  closer  to 
the  outer  wall  as  Reynolds  number  and  rad' us  ratio 
Increase.  The  same  qualitative  effect  occurs  for 
laminar  flow. 

Developing  Flow  In  the  Entrance 


Kotlow  (1985)  also  tabulates  values  of  the 
dimensionless  Inner  and  outer  wall  shear  stresses. 
The  inner  stress  Is  the  higher  and,  for 
fully-developed  flow,  these  stresses  may  be 
estimated  by  the  following  curve-fit  formulas: 


1.005  b1-5^0-2195  1O810Rev 


(12) 


Experimentally,  tne  pressure  and  shear  stress 
in  the  entrance  region  are  dependent  upon  exact 
entrance  conditions:  sharp-edged,  rounded-edge, 
bellmouth,  etc.  Some  experimental  results  for  a 
fixed  annulus  are  given  by  Oklishi  and  Serovy 
(1967),  Rothfus  et  al .  (1955),  and  Olson  and 
Sparrow  (1963).  These  results  are  discussed  by 
Kotlow  (1985). 


ffd(bn) 

9(b*8) 


(13) 


where  ff^  is  to  be  estimated  Trom  Eq.  (11). 
Equations  (12)  and  (13)  have  an  accuracy  of  about 
t91  over  the  range  Rev  *  10®-10®  and  b  -  1.01-2.Q. 


The  present  analysis  follows  Sud  and  Chaddock 
(1981)  by  assuming  that  the  entrance  velocity 
profile  ts  uniform,  u  -  Uv ,  as  in  Fig.  2. 

Boundary  layers  6(  and  i0  grow  on  the  walls  until 
they  meet  at  x  «  L* ,  the  entrance  length.  At  any 
x,  the  core  velocity  Uj(x)  is  assumed  uniform. 
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b  =  r  /r . 
o  1 


«  Figure  4.  Computed  position  of  maximum 

velocity  in  the  annulus  for 
fully-developed  flow. 


At  any  given  x,  variables  o'  and  y*  are 
scaled  with  local  shears  t^Cx)  and  t0(x)  and  the 
boundary  layer  profiles  uj  and  u0  computed  by 
Integration  of  Eqs.  (5)  and  (6).  Since  t(y)  Is  not 
known  In  the  developing  region,  the  linear  profile 
assumption  of  Wilson  and  Medwell  (1971)  was 
adopted 

t(y)  -  rw(l  -  y/5)  (15) 

for  both  the  inner  and  outer  regions. 

The  study  began  by  repeating  the  computations 
of  Sud  and  Chaddock  (1981)  for  the  special  case 
given  In  their  Figures  3  and  4.  Only  qualitative 
agreement  was  obtained,  and  the  computations 
Indicated  that  local  mass  and  momentum  balances 
were  not  accurately  satisfied.  The  discrepancy. 

If  any,  was  difficult  to  analyze,  for  Sud  and 
Chaddock  used  the  transformation  technique  of 
Wilson  and  Medwell  (1971),  resulting  In  a  double 
integral  with  extremely  complex  arguments.  In  any 
case,  the  Sud  and  Chaddock  approach  was  abandoned 
and  Instead  a  local  mass  and  momentum  balance  was 
used  to  compute  Incremental  changes  In  velocity, 
pressure,  and  shear.  Developing  flow  was  thus 
computed  for  the  ranges  Rev  -  lO^-lO?  and  b  • 
1.01-2.0. 

A  momentum  balance  over  a  short  distance  ax 
of  Inner  boundary  layer  gives  Eq.  (20)  of  Sud  and 
Chaddock  (1979): 


C(V4l)2  -  r2]  ♦  t^Ax  - 


AM.  1  1  p 

-  0  4[  J  putr  dr]  -  A[  /  p^r  dr]  (16) 

r  r 

l  1 

where  Tj  is  the  average  inner  wall  shear.  An 
exactly  similar  relation  holds  for  the  outer 
boundary  layer. 


A  moment  lm  balance  of  cr.e  entire  annular  :  ;  vw 

gt  ve3 


U3tng  values  of  (ti,t0,60)  from  the  previous  step 
and  estimating  a  new  value  for  4(,  the  velocities 
U|  and  u0  were  computed  from  Eqs.  (5,6,15)  for  tne 
next  step  x»ax.  The  Integrated  average  of  the 
velocity  profile  was  compared  with  Ub  from  Eq.  (1) 
and  t j  was  adjusted  until  mass  balance  was 
achieved.  At  the  same  time,  t0  wa3  continually 
adjusted  within  an  Inner  loop  to  meet  the 
condition  U^t  -  U^q.  The  three  momentum  balances 
were  used  to  estimate  Ax  and  S0  was  adjusted  ana 
the  procedure  repeated  iteratively  until  the 
location  of  u0  matched  that  of  uj . 

The  lnvtscld-core  Bernoulli  relation,  dp  - 
-pU-dUf,  was  omitted  in  favor  of  simultaneously 
balancing  the  inner,  outer,  and  overall  momentum 
relations.  Use  of  the  Bernoulli  relation  seemed 
to  Induce  numerical  instabilities  or  errors,  such 
that  the  pressure  gradient  or  wall  shears  fell 
below  ^helr  fully-developed  values  (Kotlow  1985). 
Several  discrepancies  of  this  type  are  seen  In  the 
tabulated  results  of  Sud  and  Chaddock  (1979),  who 
used  the  Bernoulli  relation. 

A  direct  comparison  wa3  made  with  the  special 
case  computed  by  Sud  and  Chaddock  (1979,1981):  n 
-  1.83  in.  r0  -  2.44  m,  Uv  -  61  m/s,  for  air  at 

37°C  and  0.25  atm.  The  results  are  shown  in  Fig. 

5  for  the  local  pressure  drop.  The  two  results 
are  similar,  but  the  present  computations  are 
smoother  and  perhaps  3llghtly  more  accurate. 


x  -  (m) 

Figure  5.  Comparison  with  the  special  case 

commuted  by  Sud  and  Chaddock.  (Bars 
denote  numerical  uncertainty.) 
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The  pressure  drop  and  wall  shear  streses  tn 
both  the  developing  and  fully-developed  regions 
are  smoothly  varying  functions  of  Reynolds  number 
and  radius  ratio.  They  may  be  fit  to  simple 
algebraic  correlations  which  have  an  accuracy  of 
about  five  percent.  The  friction  factor  in  the 
fully-developed  region  may  be  directly  related  by 
scaling  factors  to  the  ordinary  Moody-chart 
friction  factor  for  circular  pipe  flow. 
Development  lengths  vary  from  ten  to  forty 
hydraulic  diameters,  Increasing  slowly  with  both 
Reynolds  number  and  radius  ratio. 
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